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Polarization  characteristics  of  coastal  waters  and  their  impact  on  in¬ 
water  visibility 


Jing  Zhou1,  Alberto  Tonizzo1,  Alex  Gilerson1*,  Miehael  Twardowski2,  Derie  Gray3,  Alan 
Weidemann',  Robert  Amone3,  Barry  Gross1,  Fred  Moshary1,  Sam  Ahmed 
'Optical  Remote  Sensing  Laboratory,  the  City  College  and  the  Graduate  Center  of  CUNY,  New 

York,  NY,  10031,  United  States 

Department  of  Research,  WET  Labs,  Ine.,  165  Dean  Knauss  Dr.,  Narragansett,  Rhode  Island, 

02882,  United  States 

’Naval  Research  Laboratory,  Code  7333,  Stcnnis  Spaee  Center,  Mississippi,  39529,  United  States 


ABSTRACT 

Polarization  characteristics  of  coastal  waters  were  recently  measured  during  a  cruise  on  the  R/V  “Connecticut”  in  the 
areas  of  New  York  Harbor  -  Sandy  Hook,  NJ  region  using  a  new  Stokes  vector  instrument  developed  by  the  Optieal 
Remote  Sensing  Laboratory  at  CCNY.  The  instrument  has  three  hyperspeetral  Satlantie  radianee  sensors  each  with  a 
polarizer  positioned  in  front  of  it,  with  polarization  axes  aligned  at  0,  90  and  45°.  The  measured  degrees  of  polarization 
(DOPs)  and  normalized  radiances  as  a  function  of  angle  and  wavelength  mateh  very  well  with  simulated  ones  obtained 
with  a  Monte  Carlo  radiative  transfer  code  for  the  atmosphere-oeean  system.  In  order  to  numerically  reproduce  the 
polarized  images  for  underwater  horizontal  imaging  system  the  measured  typical  underwater  polarized  radianee  was 
used  to  estimate  the  polarized  components  of  the  background  veiling  light  and  the  blurring  effects  were  modeled  by 
point  spread  functions  obtained  from  the  measured  volume  scattering  functions  from  this  cruise  and  other  typical  oceanic 
environments.  It  is  shown  that  the  visibility  ean  be  improved  for  unpolarized  target  by  placing  a  polarizer  oriented 
orthogonally  to  the  partially  polarized  direction  of  the  veiling  light  before  camera.  The  blurring  effects  strongly  depend 
on  the  small  angle  scattering  in  the  forward  directions.  For  polarized  targets  the  Monte  Carlo  simulation  of  slab 
geometry  for  polarized  pencil  light  shows  that  the  scattering  medium  with  high  g  value  has  a  very  strong  ability  to  retain 
the  polarization  status  of  the  incident  light,  which  ean  be  utilized  to  improve  the  image  contrasts  for  targets  with  very 
different  polarized  reflection  properties. 

Keywords:  Polarization,  volume  scattering  function,  visibility. 


1.  INTRODUCTION 

Study  of  the  underwater  polarized  light  field  can  lead  to  significant  improvements  of  underwater  visibility.  Polarization- 
sensitive  vision  is  well  documented  as  serving  in  navigation  and  enhancement  of  target  detection  for  many  marine 
speeies  [1].  Knowledge  of  the  underwater  polarization  pattern  can  help  us  better  understand  the  mechanisms  involved  in 
the  sensitivity  of  animals  to  the  polarization  of  light  and  the  manners  in  which  these  animals  utilize  underwater 
polarization  Inspired  by  biological  visual  systems,  many  studies  have  applied  various  polarization  techniques  in 
underwater  imaging  systems  to  improve  the  visibility  of  selected  target  features  in  a  scattering  medium  [2-5].  Since  the 
underwater  scattered  light,  or  veiling  light ,  is  the  dominant  cause  of  image  degradation,  assessing  the  effectiveness  of 
these  polarization  techniques  under  natural  illumination  also  relies  on  the  analysis  of  the  speetral  and  geometrical 
angular  dependence  of  the  underwater  polarized  light  field.  When  compared  with  standard  radianee  intensity  data,  the 
polarized  upwelling  radianee  is  more  sensitive  to  the  intrinsic  nature  of  in-water  partieulates  such  as  partiele  size,  shape 
and  refraetive  index  [6,  7];  therefore  it  carries  more  information  on  the  background  veiling  light. 

Despite  the  importance  of  the  underwater  polarization  for  marine  applications  most  of  the  previous  studies  were  based 
on  theoretical  calculations  due  to  the  lack  of  appropriate  instrumentation  to  perform  aeeurate  and  reliable  polarization 
measurements  [6,  8].  Recently  we  developed  a  hyperspeetral  and  multiangular  polarimeter  capable  of  performing  in-situ 
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measurements  of  the  underwater  polarized  radiance  [9].  In  this  paper,  typical  patterns  of  the  underwater  polarization 
obtained  with  this  instrument  during  a  cruise  in  the  coastal  areas  of  New  York  Harbor  -  Sandy  Hook,  NJ  region  will  be 
presented,  after  a  brief  description  of  the  instrumentation.  The  correctness  of  the  experimental  data  was  confirmed  by  the 
results  of  the  Monte  Carlo  radiative  transfer  simulations  for  the  atmosphere-ocean  system.  In  order  to  determine  the 
effects  of  a  scattering  medium  on  the  quality  of  underwater  imaging,  the  point-spread  function,  obtained  from  the 
measured  VSF,  and  the  measured  polarized  pattern  are  combined  to  numerically  reproduce  the  image  of  a  standard  Air 
Force  Target  under  natural  illumination  conditions.  Results  obtained  by  using  different  VSFs  measured  in  various 
oceanic  environments  were  then  compared,  focusing  on  the  impact  on  underwater  visibility. 


2.  UNDERWATER  POLARIZATION  MEASUREMENTS 


2.1  The  polarimeter  and  the  cruise 

The  polarimeter  consists  of  three  Satlantie  Hyperspectral  radiance  sensors  with  a  polarizer  attached  in  front  of  each  one. 
The  orientations  of  the  polarizers  are:  0°  (vertical),  90°  (horizontal)  and  45°.  A  stepper  motor  was  used  to  rotate  the 
sensors  to  cover  the  full  0-180°  range  of  scattering  angles,  at  5°  steps.  A  customized  Labview  program  is  responsible  for 
the  synchronization  of  the  motor  rotation  and  the  data  stream  from  the  hyperspectral  sensors  and  for  the  displaying  the 
polanzed  spectra  in  real-time.  If  L0,  Lgo  and  L4$  are  the  radiances  recorded  by  the  three  Satlantie  sensors,  then  the  DOP  is 
given  by: 


DOP  = 


^90  )  +  (2^45  “  fy90  “  ^0  ) 


(1) 


where  L  =  Lo+Lw  is  the  total  intensity.  For  normalization  purposes,  the  downwelling  irradianee,  Ed  was  also  monitored 
with  a  Satlantie  Hyperspectral  irradianee  sensor  positioned  on  deck  of  the  boat.  Water  optical  properties  were  measured 
by  an  AC-9  instrument  (WET  Labs,  lne.).  VSF  measurements  were  obtained  with  a  custom  device  called  the  MASCOT, 
which  uses  a  658  nm  laser  diode  source  and  17  independent  detectors  to  measure  the  VSF  from  10  to  170°  in  10° 
increments.  Data  were  collected  at  8  stations  during  a  recent  cruise  on  the  R/V  “Connecticut”  in  the  coastal  areas  of  New 
York  Harbor  -  Sandy  Hook,  NJ  region,  on  July  21-23  2008. 

2.2  Typical  angular  DOP  distribution 

From  in-water  optical  measurements  at  all  stations,  chlorophyll  concentrations  were  estimated  to  be  in  the  range  1.3- 
4.8|ig/l,  minerals  concentrations  2.0-3.9mg/l,  and  CDOM  absorption  at  400  nm  approximately  0.5  m1.  Fig.  1  shows  the 
typical  angular  distributions  of  the  normalized  radiance  (Fig.  la)  and  the  DOP  (Fig.  lb)  at  510nm  for  measurements  lm 
below  the  water  surface.  Generally  speaking,  the  DOP  presents  a  bell-shaped  angular  distribution  with  maximal  values 
around  100°  and  minima  in  proximity  of  0  and  180°.  The  highest  value  of  the  DOP  is  approximately  0.4  during  sunny 
days,  which  is  lower  than  half  the  value  predicted  by  Rayleigh  theory.  This  value  can  be  even  lower  during  cloudy  days. 
Interested  readers  can  find  more  detailed  descriptions  of  the  instruments  and  data  in  [9]. 


Figure  1  also  shows  the  simulated  radiance  and  DOP  vs.  scattering  angle  obtained  with  a  Monte  Carlo  method  using  the 
measured  IOP  and  VSF  values.  As  it  can  be  seen,  simulations  and  experimental  results  match  quite  well.  Figure  2  shows 
the  comparison  between  MASCOT  and  Petzold’s  phase  function,  for  Station  1.  Similar  results  were  obtained  for  Stations 
4,  5  and  7,  due  to  similar  water  compositions. 


Staton  V  >.=5l0nm 


Scattering  Angle,  0^  (*) 


Staton  1.  X=510nm 


Fig.  1.  a)  Angular  radiance  normalized  to  downwelling  irradianee.  b)  DOP  angular  distribution  at  k  ~510nm. 
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Station  1 


Fig.  2.  Comparison  of  MASCOT  measurements  and  standard  Petzold  functions. 


3.  UNDERWATER  VISIBILITY  ASSESSMENT  FOR  UNPOLARIZED  TARGETS 


3.1  Image  modeling 

Under  natural  illumination  an  imaging  system  measures  light  from  the  target  and  additional  light  (veiling  light)  from  the 
intervening  medium  In  addition,  light  from  the  target  experiences  scattering  when  it  passes  through  the  scattering 
medium  between  the  target  and  the  camera,  hence  the  image  gets  blurred.  This  blurring  effect  can  be  fully  described  by 
the  PSF  of  the  scattering  medium  [10].  In  this  section  an  image  model  will  be  first  introduced  to  combine  together 
blurring  and  veiling  light  in  order  to  simulate  the  imaging  quality  in  the  water  scattering  medium  for  an  unpolarized 
target.  Then  the  underwater  visibility  will  be  assessed  and  discussed  by  applying  various  particulate  phase  functions  as 
well  as  various  polarization  techniques. 

Let’s  consider  an  underwater  object  illuminated  by  sun  light  and  a  camera  taking  pictures  looking  horizontally  (as  shown 
in  Fig.  3).  On  the  image  plane,  the  irradiance  signal  from  the  target  S  is  given  by  the  convolution  of  its  reflection  map 
r(x,y)  and  the  PSF  of  the  water  medium  p(x,y),  i.e.: 


S  =  r(x,y)  <8>  p(x,y) 


ty 


Camera 

*/ 


> 


z 


(2) 


Note  here  that  the  target  is  assumed  to  completely  depolarize  the  incident  light  and  diffusely  reflect  the  sun  light  acting 
like  a  lambertian  surface.  The  scattered  light  from  each  layer  dz  between  the  target  and  the  camera  reaches  the  camera 
experiencing  an  attenuation  equal  to  exp(-cz).  The  total  veiling  light  is  then  obtained  integrating  over  all  the  layers 
between  0  and  Z  (the  distance  between  the  target  and  the  camera).  Therefore  approximated  veiling  light  irradiance 
should  have  the  following  form,  assuming  that  single  scattering  is  dominant 


V (Z)  =  V {6,  oo)(l  -  exp(-cZ)) 


(3) 
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The  term  V(6,oo)  is  the  background  light  without  the  target  or  with  the  target  extended  to  infinity,  which  can  be  estimated 
from  our  angular  underwater  radiance  measurement: 

V(@,oo)  =  2?r  ^L(6)s\n0d6  ^ 

4* 

0mi„  and  0max  are,  respectively,  the  minimum  and  maximum  angle  of  the  particulate  scattering  extended  to  the  camera 
aperture.  Assuming  that  the  scattered  light  keeps  its  polarization  when  passing  through  the  medium,  the  linear  polarized 
component  of  Vom(0,oo)  has  a  form  similar  to  (4),  if  L(6)  is  replaced  by  its  corresponding  measured  linear  polarized 
components  for  0  and  90°,  which  correspond  to  horizontal  and  vertical  polarization: 


^0/90 (^»°°)  “ 


'0/90 


(#)sin  OdO 


(5) 


Then  total  signal  collected  by  the  camera  is  the  sum  of  the  blurred  image  and  the  veiling  light,  ET  =  S+V. 

3.2  Results 

A  Monte  Carlo  simulation  is  performed  by  using  the  measured  phase  function  0(0)  to  obtain  the  PSF  at  optical  depth  r  = 
5  and  single  scattering  albedo  co  =  0.9.  In  addition  to  the  measured  0(6)  from  our  cruise,  which  resembles  the  standard 
Petzold  function  (shown  in  Fig.  2),  0(0)  values  obtained  from  other  typical  ocean  environments  were  also  used  for 
comparison  purposes  Specifically,  0(6)  values  were  collected  in  Monterey  Bay,  in  waters  with  mostly  biological 
particulates,  in  Puget  Sound  with  mostly  sediment  particles  and  in  the  Hudson  Plume  which  belongs  to  oceanic 
environments  [11].  The  phase  functions  and  corresponding  PSFs  were  plotted  in  Fig.  4  together  with  the  Petzold 
functions.  The  asymmetry  factors  g ,  defined  as  the  cosine  value  weighted  by  phase  function,  are  also  calculated: 

g=  jeos  (6)P(9)dQ.  (6) 

dCl 

The  g  values  for  Petzold,  biological,  oceanic  and  sediment  cases  are  0.91,  0.97,  0.96,  and  0.94  respectively.  Although 
the  four  types  of  phase  functions  have  a  very  close  g  values  (above  0.9)  oceanic  case  is  the  least  peaked  in  the  near 
forward  direction.  They  are  angularly  flat  in  the  backward  direction  but  with  different  backscattering  levels  because  of 
the  different  refractive  indices  of  the  dominant  particulates.  For  example,  the  low  refractive  index  of  the  biological 
particles  leads  to  the  lowest  level  of  backscattering.  Typical  polarized  angular  radiances  (corresponding  to  the  angles  at 
which  the  DOP  reaches  its  maximal  value)  obtained  from  field  measurements  at  510nm  were  used  in  eq.  (5)  to  estimate 
the  veiling  light  background 


Fig.  4.  a)  Phase  functions  and  b)  PSFs  (r  -  5  and  co  =  0.9)  for  four  types  of  typical  ocean  waters. 
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When  acquiring  images  of  an  unpolarized  target  looking  in  the  horizontal  direction  (under  sun  illumination),  a  polarizer 
can  be  placed  in  front  of  the  camera  to  remove  all  or  part  of  the  veiling  light  depending  on  the  DOP  of  the  veiling  light 
itself.  This  is  due  to  the  fact  that,  in  this  geometry,  the  scattered  light  at  an  angle  corresponding  to  the  horizontal 
direction  (around  100°  scattering  angle)  is  partially  linearly  polarized,  as  it  was  shown  in  previous  section.  Fig.  5 
displays  the  results  of  the  modeled  image  for  the  four  types  of  phase  functions  shown  in  Fig.  4,  the  target  has  a 
maximum  reflection  equal  to  5%. 

The  blurred  images  under  oceanic  environment  shows  the  worst  contrast  among  the  four  types  of  water  conditions.  This 
is  related  to  the  fact  that  its  phase  function  is  the  least  forward  peaked  among  these  four  environments.  The  vertically 
polarized  images  always  show  a  better  contrast  than  the  intensity  and  the  horizontally  polarized  ones.  The  images  of  the 
DOP  are  negatives  of  the  intensity  images  because  the  more  unpolarizcd  light  from  object  represents  lower  total  DOP. 
Of  course  the  effectiveness  of  image  improvement  increases  w  ith  higher  DOP  of  the  veiling  light. 


Fig.  5.  Image  modeling  results  of  a  standard  Air  Force  Target  using  four  types  of  phase  functions  with  r  ^  5 
and  co  =  0.9). 

4.  POLARIZED  LIGHT  TRAVELLING  THROUGH  A  SCATTERING  MEDIUM 

The  utilization  of  polarized  techniques  to  enhance  underwater  detection  requires  that  the  polarized  light  produced  by 
illumination  sources  or  interactions  such  as  scattering,  surface  reflection  and  transmission  keep  its  polarization  when  it 
travels  through  a  scattering  medium.  The  visibility  improvement  based  on  polarization  is  determined  by  how  far  the 
polarization  information  can  be  carried  through  the  scattering  medium.  In  order  to  examine  how  linear  polarized  light 
spreads  and  depolarizes  when  travelling  through  the  scattering  medium,  a  Monte  Carlo  simulation  based  on  slab 
geometry  (Fig.  6)  is  performed.  The  polarized  light  is  incident  along  the  z  axis  and  the  detector  array  is  placed  in  the  xy 
plane.  The  Stokes  vector  referring  to  the  xz  plane  was  calculated  for  each  array  element. 
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Fig. 6  Slab  geomelry. 


The  phase  function  is  obtained  based  on  Mie  theory  with  Junge  size  distribution 

fir)  -  r~y  and  refractive  index  1.15  (relative  to  water)  The  following  three  cases  were  considered: 

OKin  =0.1/uw^max  =\0fjm,y  =  3,g  =  0.92 
(2  )rmin  =  0.  l/tfn,rmax  =  \0^m,y  =  5,g  =  0.82 
(3)rmin  =  0.01///w,rmax  =  =  5,g  =  0.71 

The  DOP  of  the  transmitted  light  is  plotted  in  Fig.  7  with  increasing  optical  depth  for  horizontally  polarized  incident 
light  for  the  above  three  cases.  The  depolarization  rate  depends  heavily  on  the  g  factor  of  the  phase  function  [12].  The 
smaller  the  g  factor  is,  the  faster  the  DOP  decreases  with  the  optical  depth.  It  is  shown  that  the  forward  scattering  has  a 
very  strong  ability  to  retain  its  original  polarization  status  even  for  optical  depth  higher  than  10  when  multiple  scattering 
is  dominant.  Taking  a  typieal  g  value  of  0.91  for  example,  DOP  is  over  90%  after  a  polarized  light  travelling  through  10 
optical  depths.  This  can  be  used  to  distinguish  objects  with  very  different  polarization  reflectance  properties. 

The  polarized  images  recorded  keeping  the  same  polarization  as  the  incident  light  spread  to  larger  angles  with  increasing 
optical  depths  as  shown  in  Fig.  8.  Its  eireularly  symmetrical  pattern  resembles  that  of  unpolarized  incident  light, 
implying  limited  resolution  improvement  by  using  polarized  light  illumination. 


Optical  depth  x 

Fig. 7.  DOP  as  a  function  of  optical  deplh  for  ihrcc  g  values 
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t=4  t  =  6  t-8 


Fig. 8.  The  spreading  images  of  the  horizontal  polarized  incident  light  along  z  axis  with  the  polarizer  oriented  horizontally 

before  the  detector  array 

5.  SUMMARY  AND  CONCLUSIONS 

Underwater  polarization  characteristics  were  measured  in  the  coastal  areas  of  New  York  Harbor  -  Sandy  Hook  using  a 
new  hyperspeetral  multiangular  polarimeter  recently  developed  by  the  Optical  Remote  Sensing  Laboratory  at  CCNY. 
The  angular  distribution  of  DOP  always  exhibits  a  bell  shape  with  a  maximal  value  of  0  4  around  100°  scattering  angle 
for  the  examined  waters  which  have  chlorophyll  concentrations  1.3-4.8pg/l,  minerals  concentrations  12.0-3.9mg/l,  and 
CDOM  absorption  at  400  nm  approximately  0.5  m  These  experimental  results  are  also  confirmed  with  vector  radiative 
transfer  simulations  based  on  Monte  Carlo  methods. 

The  measured  polarized  radiance  and  volume  scattering  data  were  utilized  to  model  the  underwater  images  of 
unpolanzed  targets.  The  veiling  light  level  is  estimated  from  the  measured  polarized  radiance  and  the  PSF  is  obtained 
through  a  Monte  Carlo  method  using  the  measured  VSF.  Taking  advantage  of  the  targets’  laek  of  polarization  and 
partially  horizontal  polarization  of  the  background  radiance,  the  veiling  light  can  be  reduced  by  placing  a  vertically 
oriented  polarizer  before  the  imaging  system  and  thus  improving  the  visibility.  Four  types  of  VSFs  under  typical  ocean 
environments  were  applied  to  model  the  blurred  images.  At  the  same  optical  depth  the  quality  of  the  underwater  images 
strongly  depends  on  the  small  angle  scattering.  The  ocean  condition  with  the  least  peaked  forward  scattering,  the  oceanic 
ease,  has  the  worst  image  quality. 

On  the  other  hand,  the  Monte  Carlo  simulation  of  slab  geometry  for  polarized  pencil  light  shows  that  the  scattering 
medium  with  high  £  value  has  a  very  strong  ability  to  retain  the  polarization  status  of  the  incident  light,  which  can  be 
utilized  to  improve  the  image  contrasts  for  targets  with  very  different  polarized  reflection  properties.  Evaluations  of 
imaging  the  polarized  targets  by  using  polarization  techniques  should  be  the  next  steps  of  research. 
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